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^ Abstract 



The i^* — * v T and z/„ — * v s solutions to the atmospheric neutrino problem 
are compared with SuperKamiokande data. The differences between these 
solutions due to matter effects in the Earth are calculated for the ratio of 
/x-like to e-like events and for up-down flux asymmetries. These quantities 
are chosen because they are relatively insensitive to theoretical uncertainties 
in the overall neutrino flux normalisation and detection cross-sections and 
efficiencies. A \ 2 analysis using these quantities is performed yielding 3a 
ranges which are approximately given by (0.725 — 1.0, 4 x 10~ 4 — 2 x 10 -2 eV 2 ) 
and (0.74-1.0, lx 10~ 3 -2x 10~ 2 eV 2 ) for (sin 2 29, Am 2 ) for the -> v T and 
v a v s solutions, respectively. Values of Am 2 smaller than about 2 x 10 -3 
eV 2 are disfavoured for the — > v s solution, suggesting that future long 
baseline experiments should see a positive signal if this scenario is the correct 
one. 
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Atmospheric neutrino data provides important evidence for the existence of neutrino 
oscillations Atmospheric neutrinos are produced primarily from the decays of mesons 
and muons which result from interactions between the primary cosmic ray flux and air 
molecules in the Earth's upper atmosphere. The neutrino flux consists mostly of z/ e , z/ e , 
and V jj, with the muon flavour flux expected to be roughly twice as large as the electron 
flavour flux. When these neutrinos interact with matter via the charged current they pro- 
duce the corresponding charged leptons or antileptons. Atmospheric neutrino experiments 
measure the corresponding event rates. Historically, the atmospheric neutrino anomaly has 
been defined to be the discrepancy between measured values of the ratio of /i-like to e-like 
events (~ 1.2) to the predicted ratio of about 2. The SuperKamiokande experiment has 
recently produced relatively high statistics data that considerably strengthens the case for 
an atmospheric neutrino anomaly 0. In particular, the new data provide strong evidence 
for an anomalous zenith-angle dependence for multi-GeV //-like events. Furthermore, the 
pattern of this dependence is consistent with a neutrino oscillation explanation. 

Generically, the atmospheric neutrino anomaly points for its solution towards large angle 
Vp, - v e 0, — ► v T H] or Up — > v s || oscillations, where v s denotes a sterile neutrino 0. 
However, —>■ v e is now disfavoured (though not completely ruled out) by results from the 
CHOOZ reactor-based V e disappearance experiment J7|. We therefore focus on the — > v T 
and — > v s possibilities in this paper. A major goal of future atmospheric neutrino research 
should be to discriminate between these rival solutions. To this end, the purpose of this 
paper is to compare and contrast these two most favoured oscillation modes with the latest 
SuperKamiokande data. 

Though in many respects similar, the — > v T and — > v s cases are distinguishable, 
because v T interacts via the neutral current with ordinary matter whereas v s , by definition, 
does not. Neutral current interactions with the Earth affect the evolution of the + v a 
system because of the matter effect ||. The evolution of the + v T system is, by contrast, 
identical to what it would be in vacuum. A major goal of this paper is study the magnitude of 
the matter effect. We will show that it is quite important for multi-GeV events if Am 2 < 10 -2 
eV 2 and for sub-GeV events if Am 2 < 10~ 3 eV 2 . We will also perform a \ 2 analysis of the 
two cases with respect to the SuperKamiokande data. It is interesting to note that other 
ways of discriminating between v T and v s have been suggested in the literature. Reference 
H discusses how the sensitivity of SuperKamiokande to neutral current interactions can be 
used, while Ref. discusses the importance of the matter effect for high energy neutrinos 
that produce upward-going muons. Intriguingly, the MACRO experiment |T]J sees dips in the 
upward-going muon flux at zenith angles that are qualitatively consistent with expectations 
from Ref. ]T0|. The two cases may also be distinguished in the future using long baseline 



experiments, either through the sensitivity of the detector to neutral currents, or, for higher 
energy experiments, by searching for v T appearance. 

In the water-Cerenkov SuperKamiokande experiment, neutrinos are detected via the 
charged leptons a (a = e or /x) produced from neutrino scattering off nucleons in the water 
molecules: v a N — > aX, where the identity of X will be discussed below. The total number 
N(a) of charged leptons of type a produced in either of the two scenarios considered is given 
by 
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Here d 2 F a /dEd cos £ is the differential ffux of atmospheric neutrinos of type v a of energy 
E at zenith angle £. The term is the effective number of target nucleons. The function 
d 2 a a jdqd cos ip is the differential cross section for v a N — > aX scattering, where q is the 
energy of the charged lepton and ip is the scattering angle relative to the velocity vector of 
the incident v a (the azimuthal angle having been integrated over). The function P[y a — > 
u a ;E,£) is the survival probability for a z/ Q with energy E after travelling a distance L = 

\J(R + h) 2 — R 2 sin 2 £ — cos£, where i? is the radius of the Earth and ft, ~ 15 km is the 
mean altitude at which atmospheric neutrinos are produced. Finally note that <j) is the 
azimuthal angle relative to the incident neutrino direction (see Fig.l for an illustration of all 
the relevant angles). The integration over is only non-trivial when calculating zenith angle 
binned data. In order to calculate the number of a-like events for certain energy ranges 
and within certain zenith angle bins, the integration ranges in Eq.(|l|) must be truncated 
accordingly, with the direction of the charged lepton then obtained from 



cos 9 = cos £ cos ip + sin £ cos <fi sin ip, 



(2) 



where is the zenith angle of the charged lepton. 

Since not all charged current events are used in the SuperKamiokande analysis, Eq.[I| 
must be modified to incorporate the so-called detection efficiency function. This function is 
defined by 



detection efficiency = 



Number of 1 ring charged current events 
Total number of charged current events 



(3) 



For our sub-GeV analysis we use the approximation of only including quasi-elastic scattering. 
In this approximation the detection efficiency function is set to 1 because quasi-elastic scat- 
tering always lead to 1 ring events. In the multi-GeV analysis, we use a detection efficiency 
function obtained from the SuperKamiokande collaboration. 

The survival probability P(v a — > u a ; E, £) is obtained by solving the Schrodinger equation 
for neutrino evolution including matter effects. It is given by 
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where x is the distance travelled, Am 2 the difference in squared masses, 9 the vacuum 



mixing angle and v 



fJ,,T,S 



X, 



the wave-functions of the neutrinos. The quantities A TjS (x) are 



the effective potential differences generated through the matter effect: 

AJx) = 



(5) 



and, for electrically neutral terrestrial matter [|12 
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A s (x) = ^G F N n (x) = ^G F (Y n /m n )p(x), (6) 

where Gf is the Fermi constant, N n (x) is the number density of neutrons along the path of 
the neutrino, Y n ~ 0.52 is the average number of neutrons per nucleon, m n is the nucleon 
mass and p(x) is the mass density. Our numerical calculations use the density profile of 
the Earth given in Ref. [13|]. The survival probability is given by v n {V)* v ^{V) . For 



antineutrinos the sign of A s is reversed. 

The differential flux of atmospheric neutrinos d 2 F a /dEd cos£ without geomagnetic effects 



is given in [IB], but we have used the differential flux which includes geomagnetic effects 
Wh - (For other atmospheric neutrino flux calculations, see Ref. |T8| .) 

SuperKamiokande separates its data into a sub-GeV sample and a multi-GeV sample. 
For sub-GeV events, charged leptons are dominantly produced via quasi-elastic scattering: 
u a N — > aN', where N and N' are nucleons. The state X [see Eq.(|l])] is therefore identified 



with N' for these events. We use the cross-section given in Ref. [Q. The struck nucleon 



N is either a proton in hydrogen, or a bound nucleon in oxygen. In the case of hydrogen, 
the 2-body nature of quasi-elastic scattering leads to a relation between E, q and ip from 
relativistic energy-momentum conservation. One of the integrations in Eq.fll]) is therefore 
redundant for scattering off hydrogen. For scattering off a nucleon within an oxygen nucleus, 



Fermi motion and Pauli blocking effects are incorporated via the prescription in Ref. |14] . In 



this case, there is no relation between E, q and ip because of the nuclear effects. For, multi- 



GeV events we use the inclusive cross-section for u a N — > aX given in Ref. fll5| . Although 
this cross section is not completely satisfactory for calculating absolute event rates because 
it does not incorporate low Q 2 effects such as A resonance production, it is a good enough 



approximation for calculating ratios of event rates such as up-down asymmetries |L9|] and 
N(fi)/N(e) because these quantities are relatively insensitive to details. Indeed the main 
advantage in using event rate ratios is that they are relatively insensitive to uncertainties 
in the cross-sections and the neutrino fluxes. In this way 20 — 30% uncertainties in overall 
flux normalisations and cross-sections are avoided in favour of quantities with systematic 
uncertainties of only a few percent. 

We will now define the event rate ratios used in the analysis. We first define the tradi- 
tional quantity R, where 

= (N 7n)\ ' [ ) 

V 1 * fJ,/ 1 v e) \ no—osc 

The quantities N e ^ are the numbers of e-like and /i-like events, as per Eq.(|J). The numerator 
denotes numbers obtained from Eq. ([]]), while the denominator the numbers expected with 
oscillations switched off. A class of up-down flux asymmetries for a-like events is defined by 

yr/ — C^J3 71 1 Ng n )\ osc , ^ 

no— osc 

Here N~ v denotes the number of a-like events produced in the detector with zenith angle 
cos O < — rj, while N+ n denotes the analogous quantity for cos O > r], where rj is defined to 
be positive. SuperKamiokande divides the (—1,-1-1) interval in cosO into five equal bins. 
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The central bin straddles both the upper and lower hemispheres, and is thus not useful for 
up-down asymmetry analyses. We therefore choose 77 = 0.2 in order to utilise all the data 
in the other four bins. Since z/ e 's do not oscillate in the two scenarios z/ M — > v T S which we 
consider, up-down asymmetries for e-like events, Y" e a2 are predicted to equal 1. Note that 
systematic uncertainties for up-down asymmetries are smaller than for R, because the latter 
depends on the relative flux of to u e . 

In the context of the —>■ u TyS scenarios considered here, R measures the disappearance of 
z/^'s and V^s relative to z/ e 's and TVs an d to no-oscillation expectations. The up-down asym- 
metries probe the zenith angle dependences of the neutrino fluxes relative to no-oscillation 
expectations. Both classes of quantities provide pertinent information about the pattern of 
the putative oscillations while being insensitive to the systematic uncertainties discussed 
above. 

Additional important information about the oscillation pattern is supplied by the energy 
dependences of R and the Y's. For this reason we calculate separate i?'s and Y's for the 
SuperKamiokande sub-GeV and multi-GeV samples. The sub-GeV sample is defined by the 
momentum cuts 0.1 < p e /GeV < 1.33 and 0.2 < p^/GeV < 1.5 for e-like and /z-like events, 
respectively. We also consider an alternative low-energy cut which has a lower limit of 0.5 
GeV. The alternative cut enhances the effect of oscillations because the correlation of the 
produced charged leptons with the incident neutrinos is stronger for higher energies. 



Our results for the i?'s and K's are displayed in Figs. 2-6 |20j together with the preliminary 
SuperKamiokande results 0, 

R(sub - GeV) = 0.61 ± 0.03 ± 0.05, 
R(multi - GeV) = 0.67 ± 0.06 ± 0.08, 
Y°- 2 (sub - GeV) = 0.78 ± 0.06, 
Y^ 2 (multi - GeV) = 0.49 ± 0.06. (9) 

The preliminary experimental results we use correspond to 414 live days of running. Note 
that experimental results for the alternative sub-GeV sample with p e ^ > 0.5 GeV are not 
at present available. For completeness we mention that the preliminary SuperKamiokande 
results for the e-like up-down asymmetries are 

Y°- 2 (sub - GeV) = 1.13 ± 0.08, 
Y°- 2 (multi - GeV) = 0.83 ± 0.13. (10) 

Finally note that only statistical errors are given for the up-down asymmetries since they 
should be much larger than possible systematic errors. 

These results have several significant features: (i) For the sub-GeV cases, the matter 
effects become noticeable at about Am 2 = 10~ 3 eV 2 and are really very significant at 
Am 2 = 10~ 4 eV 2 . (ii) The matter effects cut in at the higher value of about Am 2 = 10~ 2 eV 2 
for the multi-GeV cases, (iii) The up-down muon asymmetries plateau between about 10 -3 
and 10~ 2 eV 2 , with the sub-GeV plateau occuring for slightly lower values of Am 2 compared 
to the multi-GeV case. For this range of Am 2 , downward travelling neutrinos do not have 
time to oscillate, whereas upward travelling neutrinos experience averaged oscillations. The 
plateau phenomenon provides a characteristic prediction for up-down asymmetries that is 
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reasonably insensitive to Am 2 , while remaining sensitive to the mixing angle. Note also 
that the multi-GeV R flattens out in this Am 2 range, for exactly the same reason, with 
the corresponding plateau for the sub-GeV R at lower Am 2 's (except that the matter effect 
destroys the plateau for the — v s scenario), (iv) Both the sub- and multi-GeV Y's fit the 
SuperKamiokande data well for Am 2 in the range 1(T 3 - 10" 2 eV 2 . The R values also fit 
the data well in this range, but with a preference for higher Am 2 values, (v) While a fairly 
large range of mixing angle values is consistent with both of the R measurements and the 
sub-GeV Y datum, the multi-GeV Y result tends to favour maximal mixing. The multi-GeV 
R measurement also tends to favour a large mixing angle, (vi) A comparison of Figs. 4 and 
5 shows that the alternative low-energy cut increases both the up-down asymmetry effect 
and the matter effect. 

We now perform a x 2 fit to these data in order to quantify the significance of the various 
competing influences discussed above. We define the x 2 function as 

(11) 

where the sum is over the sub-GeV and multi-GeV cases, the measured SuperKamiokande 
values and errors are denoted by the superscript "SK" and the theoretical predictions for 
the quantities are labelled by "th". The r\ = 0.2 choice is understood for the up-down 
asymmetries. We include both the e-like and the /i-like up-down asymmetries in the fit. 
There are 6 pieces of data in x 2 an d 2 adjustable parameters, Am 2 and sin 2 29, leaving 4 
degrees of freedom. 

The statistical procedure we employ is approximate in the sense that ratios of Gaussian- 
distributed quantities are only approximately Gaussian themselves [21|| . The validity of 
using ratios increases as the fractional errors decrease. Since SuperKamiokande is a high 
statistics experiment, our procedure is accurate within a sufficiently small region around the 
best fit point provided this point gives a good fit (we estimate that it is approximately valid 
within the 3a region around the best fit). 

An alternative x 2 analysis can be performed by using absolute event rates rather than 
ratios. However in that case the numerical validity of the results is limited by the correctness 
of the cross-sections used. These analyses typically incorporate a 20 — 30% uncertainty in 
the event rates due to uncertain fluxes by introducing theoretical errors in addition to 
measurement errors. Unfortunately, existing analyses do not address the issue of uncertain 
cross-sections. Our analysis avoids this problem, and is therefore complementary to the 
absolute event rate type of analysis. Our work also extends other recent fits P,|2"2"|j by 
considering the — > v s case. 

The results of the x 2 fits are displayed in Figs. 7-12. Figure 7 shows the allowed region 
of (sin 29, Am 2 ) at various confidence levels for the — > v T scenario. Maximal mixing 
provides the best fit, and Am 2 values in the 10~ 3 to 10~ 2 eV 2 range are favoured. Note that 
the confidence levels are defined in the usual way by 

X 2 = Xlin + A X 2 (12) 
where A% 2 = 2.3,4.6,6.2, 11.8 for the la, 90% C.L., 2a and 3a allowed region respectively. 
Our Xmin f° r v n — > v r oscillations is Xmin = 4.5 for 4 degrees of freedom. This is quite a 
good fit to the data (allowed at the 35% level). 
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In Figure 8 we show the allowed region considering just the asymmetries instead of using 
both the asymmetries and the R ratios. This is of interest because systematic uncertainties 
for Y's are smaller than those for i?'s. Note that in this case there are 4 data points and 2 
free parameters which gives 2 degrees of freedom. 

According to Fig. 9, x 2 does not experience a deep minimum at the best fit point with 
respect to Am 2 . This reflects the plateau phenomenon discussed earlier, and shows that 
this type of atmospheric neutrino analysis will not be able to pinpoint Am 2 very precisely. 
Note that the minimum becomes shallower when the i?'s are excluded from the fit. This is 
because the plateau in R is not as pronounced as that in Y. 

Figures 10-12 show the corresponding results for the — > v s scenario. Smaller values 
of Am 2 are disfavoured in this case because the matter effect moves both R and Y away 
from the measured values. However, an order of magnitude spread in Am 2 is nonetheless 
permitted at the 3a level. It is interesting to note that present data tend to predict a positive 
signal for future long baseline experiments for the — > v s scenario, whereas the — > v T 
scenario permits Am 2 values that are too small to be probed in this manner. The value 
of Xmin f° r the — > v s scenario is xLin = 5.1 for 4 degrees of freedom. This is similar to 

— v T case and also represents quite a good fit (which is allowed at 28%). 

In conclusion, we have demonstrated that matter effects in the Earth have a significant 
role to play in comparing and contrasting the — > v T and — > v s solutions to the 
atmospheric neutrino anomaly with SuperKamiokande data. The matter effects increase 
both the ratio of //-like to e-like events, and //-type up-down asymmetries, for the — > v s 
case relative to the — > v T case for sufficiently small values of Am 2 (< 10~ 2 and < 10~ 3 
eV 2 for the multi-GeV and sub-GeV samples, respectively). Smaller Am 2 values (~ 2xl0" d 
eV 2 ) are disfavoured for the — > v s scenario, with interesting implications for future long 
baseline experiments. 
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Figure Captions 

Figure 1. The parameterisation of angles in the interaction u a N — > aX. 

Figure 2. The sub-GeV R as a function of Am 2 for various values of sin 2 29. The usual 
SuperKamiokande momentum cuts have been employed. The solid (dashed) lines pertain 
to the — > u T (z/ M — > u s ) scenario. Going from the top to the bottom curves, sin 2 29 takes 
the values 0.7, 0.8, 0.9 and 1. Note the significance of the matter effect for Am 2 < 10~ 3 
eV 2 . The dashed-dotted lines denote the preliminary SuperKamiokande result within a ±lcx 
band after 414 live days of running. 

Figure 3. The multi-GeV R as a function of Am 2 . Notation as for Fig.2. Note the signifi- 
cance of the matter effect for Am 2 < 10~ 2 eV 2 . 

Figure 4. The up-down /i-type asymmetry Y® 2 as a function of Am 2 for the sub-GeV sample 
with the usual SuperKamiokande momentum cuts. Notation as for Fig.2, except that in the 
Y > 1 region the order of the sin 2 29 values is reversed. Note the significance of the matter 
effect for Am 2 < 10~ 3 eV 2 . 

Figure 5. The up-down /i-type asymmetry Y®- 2 as a function of Am 2 for the sub-GeV 
sample with the alternative lower limit > 0.5 GeV. Notation as for Fig.2, except that 
in the Y > 1 region the order of the sin 2 29 values is reversed. Note the significance of the 
matter effect for Am 2 < 10~ 3 eV 2 . 

Figure 6. The up-down /i-type asymmetry Y®- 2 as a function of Am 2 for the multi-GeV 
sample. Both fully-contained and partially-contained events are included. Notation as for 
Fig.2. Note the significance of the matter effect for Am 2 < 10~ 2 eV 2 . 

Figure 7. The allowed region in the (sin 2 29, Am 2 ) plane for the — > u T scenario. 

Figure 8. As for Fig. 7 but with R data excluded from the fit. 

Figure 9. x 2 as a function of Am 2 along the sin 2 29 = 1 line for the — > u T scenario. Note 
the shallow minimum. Note also that the minimum becomes shallower still if R is excluded 
from the fit. The insensitivity to Am 2 corresponds to the plateau features in Figs. 2-6 (see 
text). 

Figure 10. The allowed region in the (sin 2 29, Am 2 ) plane for the — > u s scenario. 
Figure 11. As for Figure 10 but with R data excluded from the fit. 
Figure 12. As for Fig. 9 but for the i/ M — > v s scenario. 
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